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Abstract 

Corrections of two-photon exchange to proton and neutron electromagnetic form factors are em- 
ployed to study the effect of two-photon exchange on the deuteron electromagnetic form factors. 
Numerical results of the effect are given. It is suggested to test the effect in the measurement of P z 
in a small angle limit. 
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1 Introduction 

Proton and deuteron electromagnetic (EM) form factors have been studied for a long time by unpolar- 
ized electron-proton (ep) and electron-deuteron (eD) elastic scatterings and by the Rosenbluth separation 
[I] based on one-photon-exchange (OPE) . The differential cross section of the ep scattering is 

da =A Q (T N G 2 M (Q 2 )+eG 2 E (Q 2 )), (1) 

with Aq depending on kinematic variables, rjv = <5 2 /4M^ and e — [l + 2(1 + Tpf)tan 2 (9 / '2)] 1 (Mjv and 
9 are the nucleon mass and the electron scattering angle). For a long time, the extracted Q 2 -dependences 
of the nucleon EM form factors are believed to be simple dipole forms. For the proton G P E M , one 
conventionally gets 

G E (Q 2 ) = G p M (Q 2 )/f, p ~ 1/(1 + Q 2 (GeV 2 ) /0.71) 2 , (2) 

where fj, p = 2.79 is the proton magneton. Recently, the new experiments of the polarized ep elastic 
scattering were precisely carried out. The polarization transfer scattering experiments, e + p — > e + p , 
show that the ratio R p = n P G E (Q 2 )/G p M (Q 2 ) ~ 1 - 0.158Q 2 [2]. It means that R p is no longer a simple 
constant. It monotonously decreases with the increasing of Q 2 . This new phenomenon contradicts to the 



traditional knowledge shown in cq. (2). 



One way to resolve this discrepancy is to take the effect of the two-photon-exchange (TPE) into ac- 
count [31 HJ El [8] . Usually, it is believed that TPE is strongly suppressed by (Xem {oiem = 1/137). 
However, it was argued [H] that due to the very steep decreasing of the nucleon EM form factors, the 
TPE process, where the Q 2 is equally shared by the two exchanging photons, may be compatible to the 
OPE one. Some calculations of the TPE corrections to the ep elastic scattering have been done recently 
[21 IU [6l [3 [8] . There were also several other works about the TPE effect on the proton charge radius and 
on the parity-violating [TOJ, [TT] in ep scattering. The TPE corrections to the deuteron (spin 1 particle) 
EM form factors and to the e ± p processes have been also discussed in Rcfs. [T2l I13| . 

To consider the electron nucleon elastic scattering with e(k) + N(p) — > e(k') + N(p') in OPE (C = — 1), 
we have 



M? N = ^u(k')ru(k)xu{p')T^u(p), 



2 



7,<(Q 2 ) + ^^(Q 2 ) , (3) 

IAIN 



where F-j* and F% are the conventional Dirac and Pauli form factors of the nucleon. It should be men- 
tioned that in OPE, the form factors F^ 2 are real and the functions of Q 2 only. 

If TPE is considered, parity, charge- conjugation and helicity invariances lead a general expression 
Mf N = ^Wk') 1 ,u{k)xu(p^F^{Q 2 ,e)+ l ^F^{Q 2 ,e)]u{p) 



Q 



2M 



N 



+u(k'h^ 5 u(k) x w(p')7VG^(Q 2 , e)u(p) . (4) 



Here, unlike eq. (3), F^ 2 (Q 2 ,£) and G^(Q 2 ,e) are functions of Q 2 and e. F^ 2 (Q 2 ,e) can be separately 
expressed in terms of the contributions of OPE and TPE 

F lt2 (Q 2 ,e) = F h2 (Q 2 )+F$(Q 2 ,e), G A (Q 2 , e) = G { ? (Q 2 ,e). (5) 

Moreover, according to a general relation [B] 

-u(k')j ■ Pu(k) x u(p')~f ■ Ku(p) = — - — u(k')~f^u(k) x u(j/)7 M u(p) 

+ t -u{k') 1 ^ 5 u{k) x u(p')Y*7 8 «(p)> (6) 



where K = k' + k, P = p' + p, and s = (p + k) and u — (p — k') are Mandelstam variables, we have 

Q 2 

rjT = [l,F[ N (Q 2 ,e) + l -^F 2 N (Q 2 ,e) + ^F^(Q 2 ,e)f<P^ (7) 



where F$ N = with t = q 2 and 



Fi N (Q 2 ^) = h{Q 2 ,e) - °—^F 3 N (Q 2 ,e). (8) 
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Figure 1: Deuteron mass operator 



Usually a deuteron is regarded as a weekly bound system of a proton and a neutron (see Fig. 1). Many 
calculations for the EM form factors of the deuteron have been performed in different approaches in the 
literature [TU \TE[ [T|)l Ej . Recent calculations based on an effective Lagrangian approach [TH1 \T§\ have 
shown that one can reasonably explain the deuteron EM form factors with phenomenological including 
two-body operators. Note that the deuteron EM form factors receive the TPE corrections from many 
different sources in the effective Lagrangian approach. For example, the two photons directly couple to 
one of the nucleons (see Fig. 2), or the two photons directly couple to one of the two contact points A 
and B (the contact points mean the coupling points of the deuteron and its composite pn), There are also 
several other interferences between the different OPE couplings. For instant, one photon couples to A(or 
B) and another to one of the nucleons, or two photons respectively couple to the two contact points. 

In this paper, we'll study the TPE effect on the deuteron EM form factors. The TPE corrections to 
the EM form factors of the proton and neutron from the work of Blunden, Melnitchouk and Tjon [5] will 
be employed. We know that in the deuteron EM form factors the contribution from the direct coupling 
of the photon to one of the nucleons is more important than the one from the coupling of a photon to the 
contact point [T§], and the latter coupling is needed in order to guarantee gauge invariance. Therefore, it 
is expected that the TPE effect on the deuteron EM form factors is dominated by the TPE corrections 
to the EM form factors of the nucleon (see Fig. 2), and we, as the first step, consider the effect of Fig. 2 
on the deuteron EM form factors. This paper is organized as follows. In section 2 the TPE corrections 
in the eD elastic scattering are briefly discussed. Numerical results for the corrections to the EM form 
factors of the deuteron are displayed in section 3. In Sect. 4, the conclusions will be given. 

2 Two-Photon-Exchange in the eD elastic scattering 

In eD case, the electromagnetic form factors of the deuteron are defined by the matrix elements of the 
electromagnetic current Jfj,(x) according to the OPE approximation 



where p' D , A' (or pu, £, A) denote the momentum, helicity, and polarization vector of the final (or initial) 
deuteron, respectively. In eq. (9) q = p' D — pn is the photon momentum, P = pd + p'd^ Q 2 = — I 2 IS the 
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(9) 
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Figure 2: Two-photon-exchange effect on deuteron form factors from the contributions of (a) and r™ (b). 

four-momentum transfer squared, Mjj is the deuteron mass, and en is the charge of the deuteron. In 
the one-photon exchange approximation, the differential cross section of the unpolarized elastic electron- 
deuteron scattering e(fc, si) + D(j)n, £) — > e(k', S3) + D(p' D , £') in the laboratory frame is [3U] 

g?<t da 



I (OPE), Iq(OPE) = A(Q Z ) + B(Q z )tan' 



Mot! 



2' 



(10) 



where 6 is the scattering angle of the electron, (da/dfl)Mott is the Mott cross section for a structure-less 
particle with recoil effect, and the two structure functions 



A(Q 2 ) = G 2 C (Q 2 



<-t d G 2 m (Q 2 )+°-t 2 d G 2 q(Q 2 \ 



B(Q 2 ) = -t d (1 + t d )G 2 m (Q 2 



(11) 



In the above eqs. td — Q 2 /4M|,. Gj\/, Gc and Gq are the deuteron magnetic, charge and quadrupole 
form factors, respectively. They can be expressed, in terms of G\, Gi and G3, as 



Gm — G2, Gq — Gi — G2 + (1 + T£))G3, Gc — G\ + -tjjGq. 



(12) 



The normalizations of the three form factors are G c (0) = 1, Gq(0) = M^Qo — 25.83, and Gm(0) = 
1.714. Note that in eqs. (10) and (12), there are two unpolarized structure functions A and B, and 
three independent form factors Gc, Gq and Gm for the deuteron. To determine the three form factors 
completely, one needs, at least, one polarization observable. The optimal choice is the polarization T20 
(or P zz ) [H]. 



Considering both OPE (C = —1) and TPE (C = +1), and taking the Lorentz, party, and charge- 
conjugation invariances into account, one obtains the most general form of the eD elastic scattering 

mum, 
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where 
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(14) 



(15) 
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(16) 



General speaking, the form factors G[ with i = 1,6, are complex functions of s = (po 
Q 2 = —(k — k') 2 . They can be expressed as 

G' i (s,Q 2 ) = G i (Q 2 ) + G ( i 2 \s,Q 2 ), 

(2) 

where G, represents the contribution arising from the one-photon exchange, and G\ stands for the rest 
which would come mostly from TPE. In the OPE approximation, G' 4 — G' b — G' e = 0. It is easy to see 
that Gi (i = 1, 2, 3) are of order of (qem)" and G\ 2 ' (i — 1, ...6) are of order olem- Moreover, G- = G^ 
for i = 4, 5 and 6. 



To take the TPE corrections to the proton and neutron (see Fig. 2) EM form factors into account, 
and to study the TPE effect on the EM form factors of the deuteron, we directly calculate the matrix 
element of < p' D , A' | J£(0) + J™(0) \p D ,X>, where 

J£<"(0) =| pn X pn | f>>"(0) \pn><pn\. (17) 

The effective interaction between the deuteron and its composites (pn) is [19] 

C D = 9dD^+(x) J dy^ D (y 2 )p(x + \y)C lfl n(x - + H.c, (18) 

where G is the charge conjugate matrix. The correlation function characterizes the finite size of the 
deuteron as a pn bound state and depends on the relative Jacobi coordinate y, in addition, x being the 
center-of-mass (CM) coordinate. The Fourier transformation of the correlation function reads 

<Mj/ 2 ) = /(04 e_w ^(V). (19) 

A basic requirement for the choice of an explicit form of the correlation function is that it vanishes suffi- 
ciently fast in the ultraviolet region of Euclidean space to render the Feynman diagrams ultraviolet finite. 
We adopt a Gaussian form, §d{Pe) = ex P( — Pe/^d) > f° r the vertex function, where pe is the Euclidean 
Jacobi momentum. Here, Ad is a size parameter, which characterizes the distribution of the constituents 
inside the deuteron. 

The coupling go of < pu, A | pn >= *(A) is determined by the compositeness condition [231 12"4"1 
1251 126j . The condition implies that the renormalization constant of the hadron wave function is set equal 
to zero: 

Z D = 1-E' D (M D ) = 0. (20) 
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Here, S^,(M|,) = g 2 y Il' D (M^ ) ) is the derivative of the transverse part of the mass operator X^f, which 
conventionally splits into the transverse £_d and longitudinal parts as: 

E^(p) - fff S D (p 2 ) + ^^(p 2 ) , (21) 

where 

gf = g^-p^/p\ gfp a = Q. (22) 

The mass operator of the deuteron in our approach is described by Fig. 1. If the size parameter hp is 
fixed, the coupling gu is then obtained according to the compositeness condition eq. (20). 

Note that the current of photon-nucleon with TPE has an additional structure F 3 N as shown in eq. 
(7). An explicit calculation of the matrix element including TPE (see Fig. 2) gives 

Mf D = ^S(fc0 7 ^(fc)C(AO[^ + «]ep(^)' (23) 

where 

xTr (j°S F (k + pf D )f*S F (k + PD )-fS F (k)j (24) 
is the contribution from Fig. 2(a) and 

xTr(~t°S F (k)YS F {k - PD )r™S F (k - pf D )j (25) 

is the one from Fig. 2(b). From the explicit expressions of eqs. (24-25), it is found that the form factor 
of F^ 2 contributes to the charge, magnetic and quadrupole form factors of the deuteron. When the TPE 
effect is included, it provides a new form factor of the nucleon, F 3 N . This new form factor contributes, 
in our approach, to the charge and quadrupole form factors of the deuteron. Particularly it also gives a 
contribution to G' 6 — G^ with the structure of Mq and q ■ M 6 = 0. The explicit expression for the new 
form factor of G' 6 contributed by F 3 P is 

n'P 2 f°° dadfidj —prl^+M-'p 2 p 2 

G e ~ 9 ^/ ( 47r )2 A 2^ 2 e F/(Q ,e)H^(a,(3,j,Q ,e), (26) 



where 



with 



M 2 

+4— 2- {-2 Xl xl +xi- x\x 2 + x\- x\- x z 2 ) (27) 

N 



2 + « + /3 + 7 , A = ^ 



6 



and 



d = (a + f3 + ~f)M 2 N -(^ + a + f3)M 2 D 

d 2 = (l + a + /3) 2 Af2+(i + a )(i+/3)Q 2 , (28) 



1 1 f=° f dad/3 _{ n ,. , os.2 a + {3 + 2a(3 2 



exp _2(a + /3)/i^ r + - 7 -f (29) 



9 2 d ^ 2 Jo Jo (l + « + /3)3 *\ 2(l + a + /3) 

a „ x/ 2 1 (l + 2a)(l + 2/9) 2 

x[(a + + 2a/?) fj, 2 N + — — — + — m ; > d 



2 1 (l + 2a)(l + 2/?) 2 \ (l + 2a)(l + 2/3) \ 

2(l + a + /3) + 4(1 + a + /?) 2 to J + 2(l + a + /3) 

where /j# = ma /Ad with H — N, D. In eqs. (26) and (29), the integration variables a, /? and 7 are the 
Feynman parametrizations. The contribution from F 3 n to the new deuteron structure function G'§ can 
be obtained from eq. (25) in the same way. 



3 Numerical results 

In [27], the discrepancies between OPE and TPE have been carefully analyzed for the ep elastic scat- 
tering. The empirical estimates of the TPE amplitudes in the ep elastic scattering are given based on 
the assumptions about the angular dependence of the amplitudes which is limited by the precision of the 
Rosenbluth data [28] . It is also assumed that the entire form factor discrepancies are because of the new 
form factor F%. In a recent paper about a global analysis of the proton elastic form factor data with the 
two- photon exchange corrections [35] , the input TPE corrections are following the formalism of Blunden 
et al. [5], rather than the phenomenological corrections extracted in [57]. It is found that the value of 
Y27 is much smaller than that extracted in the phenomenological analyses. 

In our numerical estimates of the TPE effect on the EM form factors of the deuteron, we will also 
use the TPE corrections to the proton and neutron EM form factors following the formalism of Ref. 
[5]. It should be mentioned that the TPE corrections are ^-dependent (or e-dependent). One can get 
G^em/G'e m f° r the proton and neutron as well as F 3 ' n (Q 2 1 e) from the obtained Y^f 1 with 

v p { m(Q 2 ,e) \ K° ( P 3 (Q 2 ,e) \ 

Y * = Re \^w^\ J = ^r N Re {c^m ) ' (30) 

where / = M^^J±±I^/t n {1 + t n ) = \M N K°. To explicitly show the TPE effect on the deuteron EM 

(2) (2) 
form factors, we display, in Figs. 3-5, the ratios B4 = G\ /G^ with i = C,M,Q, where G\ are the 

contributions from TPE of Fig. 2, and G[ are taken from the phenomenological parametrization of the 
deuteron EM form factors [3U] as empirical data. In the three figures, we respectively choose the scat- 
tering angle 9 = 7r/6, tt/2 and 5ir/6. According to the constraint condition that the deuteron is bound 

as <| r~ 2 |>< 0.02GeV 2 [14], we select a typical parameter = 0.30 GeV. Moreover, in Figs. 6 and 
- (2) 

7, we present our predictions for the new form factor G' 6 = G 6 where the contributions from Fig. 2(a) 
and Fig. 2(b) are given, respectively. 
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Q 2 (GeV 2 ) 

Figure 5: Rq(Q 2 ,6) for three 9. 

Note that the measured EM form factors of the deuteron appear differently from those of the nucleons, 
since they have crossing points. The data indicate that G' c or G' M respectively has a crossing point at 
Q 2 ~ 0.7 GeV 2 or 2.0 GeV 2 [SU]. Fig. 5 shows that the TPE effect on the deuteron quadrupole form 
factor are greatly reduced comparing to the corrections to the proton and neutron. In Figs. 3 and 4, the 
peaks result from the crossing points of G' c and G' M . The two figures mean that the typical magnitudes 
of the ratios Rc and Rm, due to the TPE corrections in our approach, are always less than 10%, and 
the corrections to G' M are more remarkable than those to G' c . Moreover, Figs. 3-7 tell that the TPE 
corrections to G' c M q and to G 6 are ^-dependent. We also find the two contributions from the proton 
Fig. 2(a) and neutron Fig. 2(b) always cancel each other and make the total magnitude of G 6 being 
smaller. 

In fact, according to the analyses of Refs. [T21Q35], it is expected that the measurements of the single 
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polarization observables P x (Tn) and P z (Tio) are useful to test the TPE effect. Here we know that the 
contribution of G 6 to the polarization P x , is 



P.P ~ - 3 j^td Vr D (l + r D )tan(-)G M Re(G' 6 ). (31) 



Comparing to the P x in the OPE approximation 



P x = - 3 a/td(1 + T D)tan(-^)G m{Gc + ^t d Gq), (32) 



we find 



M^i Gc + 3^Gq 

The effect of G 6 on the polarization of P x is shown in Fig. 8 in the three cases of = 7r/6, 7r/2 and 57r/6. 
One sees the ratios are less than 1% in the range of 0.5 < Q 2 < 3 GeV 2 . In Fig. 8 the maximum points 
are expected to result from the minimum point of Gc + \tGq. 

Moreover, the polarization P z in OPE is 

Pz = ^v/^(l + ^)W^)G 2 M (34) 
and the contribution to P z from G 6 is 

P? ] - -| Wtd(1 + r D )G M Re(G' 6 ) (35) 

which is also ^-dependent since Gg is. In the small angle limit, the contribution of OPE to P z vanishes and 
the one from G 6 remains no-vanishing to the contrary. Therefore, it is expected that the measurement 
of this polarization, P z , in the small 9 limit can easily show the TPE effect. For the ratio, we get 

R(P \- P * 2) - It M ° Re{G ' e) ftffl 
P z tan l {\)GM 



In Fig. 9, we display the ratio of R(P Z ) in the three cases of 9. A larger TPE effect on P z than on P x is 

> 2 ' 



seen since the denominator of the ratio in eq. (36) is proportional to ian 2 (|). 



4 Conclusions 

To summarize, we have explicitly given the TPE corrections to the conventional form factors of the 
deuteron G' c M q. In our approach, the TPE corrections to the nucleoli EM form factors (see Fig. 2) are 
considered. We find that the new from factor of the nucleon with TPE, F^ p ' n , not only contributes to the 
form factors G' c q of the deuteron, but also provides a new form factor of the deuteron Gg . According 
to the formalism of Ref . [S] , we numerically estimate the TPE effect on the deuteron EM form factors, 
and we get the ^-dependences for all the TPE corrections. It is suggested the TPE effect can be tested 
in the measurement of the single polarization of P z (T±q) in the small angle limit. In addition, we find 
that the TPE corrections to G' M are more important than those to G' c q . 
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Of course, the TPE effect, we considered in this work, only results from the sources of the direct 
couplings of the two photons to one of the two nucleons inside the deuteron (Fig. 2). There are several 
other sources of TPE which could be included in our future calculation as next step. An overall estimate 
of all the TPE corrections to the deuteron form factors is in progress. 
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